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Conformal methods for massless Feynman 
integrals and large Nf methods 

J. A. Gracey 



Abstract We review the large N method of calculating high order information on the 
renormalization group functions in a quantum field theory which is based on con- 
formal integration methods. As an example these techniques are applied to a typical 
graph contributing to the j3-function of 0(N) (j) 4 theory at 0(1 /N 2 ). The possible 
future directions for the large N methods are discussed in light of the development 
of more recent techniques such as the Laporta algorithm. 



1 Introduction 

One of the main problems in renormalization theory is the construction of the renor- 
malization group functions. These govern how the parameters of a quantum field 
theory, such as the coupling constant, depend on scale. In situations where one has 
to compare with precision data, this ordinarily requires knowing the renormaliza- 
tion group functions to very high orders in a perturbative expansion. The quantum 
field theories we have in mind are not only the gauge theories of particle physics 
but also the scalar and fermionic ones which arise in condensed matter problems. 
These are central in understanding phase transitions. To attain such precision in per- 
turbative expansions means that large numbers of Feynman diagrams have to be 
determined with the number of graphs increasing with the loop order. Moreover, as 
the order increases the underlying integrals require more sophisticated methods in 
order to deduce their value analytically. The widely established methods of com- 
puting Feynman graphs will be reported elsewhere in this volume. Here we review 
an alternative approach which complements explicit perturbative techniques. It does 
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so in such a way that for low loop orders there is overlap but at orders beyond that 
already known part of the perturbative series can be deduced at all orders within a 
certain approximation. This is known as the large N or large Nf method where N is 
a parameter deriving from a symmetry of the theory such as a Lie group or the num- 
ber of massless quark flavours, Nf, in Quantum Chromodynamics (QCD). In this 
method the Feynman graphs are related to those of perturbation theory but because 
of the nature of the expansion parameter, the powers of the propagators appearing 
in such graphs are not the canonical value of unity but instead differ from unity by 
0(e) where £ corresponds to the regularizing parameter of dimensional regulariza- 
tion. In addition beyond leading order in the 1 /N expansion, the propagator powers 
will include the anomalous dimensions in addition to the leading or canonical di- 
mension. Therefore, standard perturbative techniques such as integration by parts 
requires care in its use since one may not be able to actually reduce a graph to a 
simpler topology. Instead a different technique has had to be refined and developed. 
It is based on a conformal property of Feynman integrals and we review it here in the 
context of the large N methods. Though it has had some applications in perturbative 
computations. 

The article is organized as follows. We devote the next section to the nota- 
tion and techniques of computing Feynman graphs using conformal methods in 
^-dimensions. We focus on the general two loop self energy graph in the subse- 
quent section and review the work of 0][2], upon which this review is mostly based, 
and others in the methods of evaluating it. These techniques are then applied to a 
problem in scalar quantum field theory in Section 4 where a graph with 10 internal 
integrations is evaluated exactly in cZ-dimensions. We conclude in Section 5 with 
thoughts on the direction in which the technique could be developed next given 
recent advances in the computation of Feynman graphs using conventional pertur- 
bative techniques. 



2 Notation and Elementary Techniques 

We begin by introducing the notation we will use which will be based on lUJ [2j. 
There Feynman graphs were represented in coordinate or configuration space nota- 
tion. By this we mean that in writing a Feynman integral graphically the integration 
variables are represented as the vertices. By contrast in momentum space represen- 
tation the integration variables correspond to the momenta circulating around a loop. 
So in coordinate space representation propagators are denoted by lines between two 
fixed points, as illustrated in Figure 1 . 
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Figure 1 . Coordinate space propagator. 
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There the power of the propagator is denoted by a number or symbol beside the line. 
One can map between coordinate and momentum space representation by using a 
Fourier transform. In the notation of JT] |2] we have 

1 - a[a) f ^J^ (1) 



(x 2 ) a 2 2a n^ Jk {k 2 y- a 

where x is in coordinate space and k is the conjugate momentum. Also for shorthand 
we set 

d = 2\i (2) 

which is used throughout to avoid the appearance of d/2 in the Euler F-function. 
This symbol should not be confused with the mass scale appearing in renormaliza- 
tion group equations. Clearly 

<■> - %r 

which is singular when a = ji + n where n is zero or a positive integer. Also a(a) 
vanishes at the negative integers. The elementary identity 

a{a)a{ji-a) = 1 (4) 

follows trivially as does 

( \ «(a-l) 
fl (°0 = 7 TV? — — T ( 5 ) 

from the F-function identity F(z+ 1) = zT(z). With this notation the elementary 
one loop self energy graph in momentum space is replaced by chain integration in 
coordinate space representation. This is represented graphically in Figure 2 



• • = v(a,j3,2 J u-a-j3) • 







where, (HQ, 



Figure 2. Chain integration. 



v(a,/3,y) = K^a{a)a{P)a{y) . (6) 



However, in practice Feynman graphs have more complicated integration points. In 
other words in coordinate space representation one has more than two lines inter- 
secting at a point. Therefore, more involved integration techniques are required to 
evaluate the Feynman integrals. One very useful technique is that of uniqueness or 
conformal integration which was introduced in three dimensions in (3). It has been 
developed in several ways subsequently and specifically to <i-dimensions. For ex- 
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ample, see (4|- We follow [Q~|[2) and use the rule represented in Figure 3, where z is 
the integration variable, 




Figure 3. Conformal integration when a + + y = 2/i. 

which follows when the sum of the exponents of the lines intersecting at the 3-point 
vertex add to the spacetime dimension 

a + J3 + J = 2m • (7) 

This is known as the uniqueness condition. By the same token if a graph contains a 
triangle where the lines comprising the triangle sum to fi such as 

(M-a) + (m-J3) + (/I -y) = n (8) 

as is the case in Figure 3, then the unique triangle can be replaced by the vertex 
on the left side. There are several methods to establish the uniqueness integration 
rule. If one uses standard text book methods such as Feynman parameters then the 
integral over z can be written as 

n»r(ii-a)r(a + p + Y-ii) 



r{p)r(r)rOi) 



o [b(\-b)(x-y) 2 ] a+ P + y-V 

prior to using, (Q. When that condition is set then the hypergeometric function col- 
lapses to the geometric series and allows the integration over the Feynman parameter 
b to proceed which results in 

ftT(/i-a)r(/i-j3)r(/i-y) 

r(2n-p-Y)r(p)r(Y)(y 2 n(x- y ) 2 }P+r-» 

x 2 Fi(a,ii-y-,2li-P-r,l-^) (10) 
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Applying the uniqueness condition a second time produces the right hand side of 
Figure 3 since the hypergeometric function again reduces to the geometric series. 
This is in such a way that the canonical propagators emerge. 

An alternative method is to apply a conformal transformation on the coordinates 
of the integral, @. In this approach, which is applicable to any graph in general, 
one external point is labelled as an origin and given as a coordinate. The other 
points are denoted by coordinates x, y and z. The conformal transformation changes 
the integration coordinate as well as the external points through 



Xn 

->-£. (ID 
x L 



Thus for two coordinates y and z undergoing such a transformation we have the 
lemma 

(y_ Z )2 _> £z£. (12) 

y z z z 

An integration measure also produces contributions to the lines joining to the origin 
since 

Therefore, for the vertex on the left side of Figure 3 this transformation produces 
the intermediate integral of Figure 4. 




Figure 4. Vertex of Figure 3 after a conformal transformation with base at where 

a = 2\i — a — $ — y. 

To complete the integration requires setting the uniqueness condition <(7j which pro- 
duces a chain integral since the line from to z is absent from the graph. To com- 
plete the derivation one undoes the original conformal transformations to produce 
the right hand side of Figure 3. If one compares the two derivations, the latter is in 
fact of more practical use. This is because it avoids the use of writing the original 
integral in terms of Feynman parameters which would become tedious for higher 
order cases. Also it is simple to implement graphically. 

Having recalled the derivation of the uniqueness rule it is straightforward to see 
that there is a natural extension. In the first derivation there was not a unique way to 
collapse the hypergeometric function to an elementary type of propagator. Instead 



6 



J.A. Gracey 



this will happen if the sum of the exponents is (2/x + n) where n is a positive integer. 
Athough the collapse in this case will not be to the geometric series, it will reduce to 
simple algebraic functions which are of the propagator type. So, for instance, when 
n = 1 we have the result of Figure 4, 0, 




v(a-i,p-i,y) o 

(0-1)0-1) V- P 



v(g-l,]3,y-l) 

(g-i)(y-i) 



v(g.ft-l.y-l) 
(jB-l)(y-l) 



J" 




/i - a 

Figure 5. Conformal integration when a + J3 + y = 2\l + 1. 

A similar rule has been constructed and used in [|4J. We will use Figure 4 later in 
order to simplify various integrals. 



3 Two Loop Self Energy Graph 

We can illustrate some of the techniques of conformal integration by considering the 
massless two loop self energy graph with arbitrary powers, a, on the propagators. It 
is illustrated in Figure 6 




Figure 6. Two loop self energy graph in coordinate space representation. 
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where we have used the coordinate space representation. Thus the vertices are inte- 
grated over rather than the loop momenta. To clarify, the integral of Figure 6 is 

7(0=!, 0= 2 , 0=3, 0:4,0:5) = £ (3,2)0, (( x -y)2)<*i(( x - z )2)a 3 (z 2)« 4 {{y _ z) 2)« 5 < 14 > 

where J y = J j^ji ■ The structure of this integral has been widely studied and we 
briefly highlight several properties of relevance. The analysis of (6j [TJ determined 
that the symmetry group of the graph was Z2 x $6 which has 1440 elements. Ex- 
ploiting this the e expansion of the integral in d = 4 — 2e with propagator powers of 
order e from unity was determined up to 0(e 6 ), flUE). At 0(e 5 ) it was discovered 
that the first multi-zeta value occurred, Q. Specifically 

7(1,1,^-1,1,^-1) =6C 3 +9C4£ + 7C 5 £ 5 

+ 5 -[Ce-2Ci] e 3 -i[9lC7 + 120C3C4]e 4 

+ [653440C 5 C3 - 7059417C8 + 576F 53 ] e 5 

+ 0(e 6 ) (15) 

where is the Riemann zeta function and F53 = L»>m>o tt in the original no- 
tation of 0. Subsequent to this it has been shown that the only numbers which 
appear in the full series expansion in e are mutiple zeta values, (9). While the work 
of (6]|3 illustrated the power of group theory to evaluate master integrals explicitly, 
using conformal integration allows one to relate two loop self energy integrals by 
exploiting the masslessness of the original diagrams. This was originally developed 
in IHEl an d we summarize that here as there appears to be scope nowadays to take 
this method to three and higher loop order graphs. 

The transformations developed in (fZ| fall into several classes. The first is that 
derived from the elementary use of the Fourier transform. Writing 

r 

7(0=1, 0=2, 0=3,0=4, 0:5) = ( x 2}D-2v ^ 

where D = °k an d T is independent of x and corresponds to the value of the 
integral, then taking the Fourier transform produces an integral which is also the 
two loop self energy. Though the propagator powers are different. In this sense one 
can say that the graph is self-dual which is not a property all Feynman graphs have. 
Thus, ED, 

7(0=1, 0=2, 0=3, 0=4, 0=5) = !Z} Z ' : /(M-»2,M-»3,M-«4>M-"l>M-«5) ■ (17) 

a(p — 2/i) 

This transformation is known as the momentum representation or MR. It can be 
easily generalized to other topologies and there is a simple graphical rule for this. 
Although not immediately apparent from the self energy because of the self-duality, 
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each 3-valent vertex of the original graph has an associated triangle in the dual 
graph. For other topologies 4-valent vertices are mapped to squares and 5-valent 
vertices to pentagons with a clear generalization pattern. 

A set of less obvious transformations can be deduced from the uniqueness con- 
dition. First, we define the shorthand notation, J2], 

S\ = <X\ + O1 + O5 , S2 = O3 + O4 + O5 

h = 01 + 04 + 05 , t2 = 052 + 03 + 05 (18) 

and illustrate the technique for one case. If one considers the central propagator it 
can be replaced by a chain integral. Although there are an infinite number of ways 
of doing this one can choose the exponents of the chain so that the top vertex is 
unique. In other words 

1 



((y-z) 2 )^ 

L 



v(2fi - O] - a 2 ,si - n,n - as) J u {(y-u) 2 ) 2 ^-^-{{u-z) 2 ) s '-^ 



(19) 



where u is the intermediate integration point. As the v vertex of Figure 6 is now 
unique the conformal integration rule can be used to rewrite the integral. This results 
in, E), 

. . 0(01)0(02)0(05) , . 

/(ai, 02,03,04,05) = -i — — — -I(jx -a 2 ,n-a u a 3 , o 4 ,si -fi)- (20) 

In the notation of |2) this transformation is known as f. It is elementary to see 
that there are five other such transformations which are denoted by /*, \, I, \ 
and ) /. The syntax is that when an arrow points in a general upwards direction 
it is a transformation on the y vertex and by contrast in a downwards direction it 
relates to the z vertex. The propagator which one replaces by a chain to make the 
vertex unique is in correspondence with the direction of the arrow. While these six 
transformations operate on the internal vertices there are two which act on each of 
the external vertices. One can complete the uniqueness of one of these by realizing 
that the integral itself is a propagator with power (D — 2jli) as indicated in ( TToT ). 13. 
For example, if the right external point is chosen as the base integration vertex then 
the appending propagator has power (2jtt — 02 — O3). This produces 

/ x 0(09)0(03) , . 

/(01, 02,03,04,05) = — \~ — / (ai,/i-a 3 , Ai-a 2 , 04,f 2 -ju) (21) 

a{D — 2n)a(2ix — 1\) 

and this is denoted by <—. The corresponding transformation on the left external 
point is called — K 

The final set of transformations are based on the conformal transformations (fTTT i 
and ( fT2b together with the effect they have on the two vertex measures, |2 |. One 
can choose either of the external vertices as the origin of the transformation. Once 
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decided the result of the conformal transformation is that all propagators joining to 
the origin have their powers changed to the difference of 2ji and the sum of the 
exponents at the point at the other end of that propagator. This means all points 
including those not directly connected to the base point in the first place. For the 
two loop self energy there are no such points but for higher loop graphs this will 
be the case. We will give an example of this in Section 4. Thus the conformal left 
transformation is, 

I(a l ,a2 1 a 3 ,a 4l a 5 ) = l{2\i-si,aq., a 3 ,2/i-5 2 , as) (22) 

where there is no /"-function factor and this is denoted by CL in contrast to CR 
which is the transformation based on the right external vertex as the origin of the 
conformal transformation. The full set of transformations and the result of applying 
each to the graph of Figure 6 are summarized in a Table in 0. However, as brief 
examples of the transformations the integral of (fT5t is related as follows 

7(l,l,jU — 1,1 ;J U — 1) = 7(>-l,jU-l,l,;U-l,l) 

(23) 

Though the latter follows from a simple rotation of the integral as well. 

Aside from the transformations there are other techniques which allow one to 
evaluate the two loop self energy and higher order graphs. Perhaps the most ex- 
ploited is that of integration by parts which was introduced for ( fTol in ifTOl . It de- 
termined that the first term in the e expansion of /( 1 , 1 , 1 , 1 , 1) was 6^3 and has also 
been used in other applications, J2). Indeed more recently the technique has been 
developed by Laporta in ffTTl to produce an algorithm which relates all integrals in 
a Feynman graph to a base set of master integrals. These can then be evaluated by 
direct methods to complete the overall computation. In the coordinate space repre- 
sentation we use here the basic rule is given in Figure 7 




Figure 7. Integration by parts in coordinate space representation. 
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where the + or — on a line indicates that the power of that propagator is increased 
or decreased by unity. For example, with this, ifTO) , 



7(1,1,1,1,1) 



v(l,l,2/x-2) 
{pi -2) 



[v(l,2,2M-3)-v(3-|i,2,3M-5)] (24) 



which can be expanded in powers of £. Clearly the series can only involve rationals 
and £„. Indeed the rule can also be applied to more general cases. In (2) it was shown 
that 



I(a,fi- 1,/x — l,j3,/i — 1) = 



a(2ju-2) 



a(a)a(2 - a) 
(l-0)(a + /3-2) 

a(a+j3-l)a(3-a-j3) 
+ (a-l)(/3-l) 
a(/J)a(2-/J) 



(l-o)(o + p-2) 



(25) 



for arbitrary a and /J. However, not all graphs can be integrated by parts. An exam- 
ple of such a case is 1(1 , a, J3 , y, 1 ) for non-unit a, j3 and y. Another example is (fTBI l. 
0, whose expansion has a non-Riemann zeta value at some point in the expansion. 
Indeed this is perhaps an indication of an obstruction to integrability. 

While integration by parts allows one to reduce the powers of various propagators 




»i(M-»i-l) 
(a 2 -l)(a 5 -l) 




(a 2 +ai -fi-l) 
(«2-l) 




Figure 8. Reduction formula for two loop self energy based on the generalized /• 

transformation. 
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by unity within a Feynman diagram it is not the only method to achieve this. A 
modification of the uniqueness method can be used to derive rules similar to Figure 
7. Specifically if one chooses the exponents of the propagators to be (2ju + 1) then 
one finds the extension given in Figure 5. Using this rule and repeating the analysis 
of the transformations on the two loop self energy graph provides relations specific 
to this topology, 0. For instance, extending to have the upper vertex exponents 
summing to (2/i + 1 ) gives the relation in Figure 8 where the + or — on the right side 
indicates that the exponent of that line is increased or decreased by unity. In Figure 
8 provided a.2 ^ 1 and a$ ^ \ then the powers of the respective propagators can 
be reduced by unity. However, this restriction is a drawback if one wishes to reduce 
graphs which have unit exponents. Instead it is possible to extend the method which 
produced the relation of Figure 8. For instance, rather than begin with the general 
two loop self energy and applying the generalized uniqueness rule, one can use 
one of the transformations of [2] and then apply a rule like that of Figure 8 before 
applying the transformation inverse to the original one. In this way one can build up 
a suite of relations. 




Figure 9. Another reduction formula for two loop self energy. 

One such useful relation is illustrated in Figure 9 which is derived in several 
stages. The first is to construct a relation similar to that of Figure 9 by first applying 

to the graph of Figure 6 and then undoing it by applying the rule of Figure 5 
to the same external vertex. This produces a relation where ?2 increases by unity in 
each of the three resulting graphs. The second stage is to apply this rule to the graph 
of Figure 6 after a CR transformation has been enacted. To complete the derivation 
the final step is to undo with another CR transformation. Thus the value of each 
graph on the right hand side of Figure 9 is one less than that of the graph on the left 
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side. This reduction has coefficients on the right hand side which are non-singular 
for unit propagators. Other rules can be derived by this method and a fuller set are 
recorded in Appendix B of |[T2l . It is worth noting that similar rules based on the 
generalized uniqueness where developed in J4]. 



4 QFT Application 

Having discussed the general techniques for determining massless Feynman inte- 
grals using conformal methods, we illustrate their usefulness in a practical problem 
in a quantum field theory. Specifically we focus on the determination of the criti- 
cal exponents at a phase transition in various models in the large N expansion. The 
background which we describe here is based on a series of articles, 0]|2][T3), where 
exponents were determined in li-dimensions at 0(1 /N 2 ) and 0(1 /N^). The fact that 
cZ-dimensional results are computable means that information on the renormaliza- 
tion group functions can be deduced in various spacetime dimensions. This is due 
to a special feature of critical point field theories and that is that at a non-trivial 
fixed point of the renormalization group flow the critical exponents correspond to 
the associated renormalization group function at that fixed point. Thus information 
on the renormalization group functions is encoded in these exponents. Moreover, 
at a fixed point several quantum field theories can lie in the same universality class 
despite having different structures. This is invariably as a consequence of a common 
interaction in the Lagrangian. Thus the same exponents can be used to access the 
structure of the renormalization group functions of two different theories. Further, 
as the spacetime dimension d is not used as a regulator, information on the expo- 
nents can be deduced simultaneously in several different dimensions such as three 
and four. For more background to the use of the renormalization group equation at 
near criticality in quantum field theories see, for example, lfl4l . 

For the application of the conformal methods we consider here we concentrate 
on the 0(N) nonlinear a model which is critically equivalent in ^-dimensions to 
0(N) 4 theory. For the latter theory the Lagrangian is 

L = \(d^'f + f (f'f) 2 (26) 

where g is the coupling constant and 1 < i < N. Introducing an auxiliary field a 
equates this Lagrangian to 

L = \(d^'f + ic7(ff) - ^ . (27) 
2^2 2g 

At criticality it is the interaction which drives the dynamics and thus it is straight- 
forward to see that in this formulation the Lagrangian interaction is the same as that 
of the 0(N) nonlinear a model when the fields are constrained to lie on an (N— 1)- 
dimensional sphere. The constraint would have a final term linear in a rather than a 
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quadratic one together with a different coupling constant. This essentially is the ori- 
gin of both field theories being in the same universality class. The linear or quadratic 
terms in a at criticality serve effectively to define the structure of the propagators. 
In coordinate space representation these are, |Q]|2], 

sij 'a d 

(*W(*)> = ^55 > <«(0M*)> = pfi (28) 

where A and B are x-independent amplitudes and a and j3 are the scaling dimensions 
of the fields. The latter comprise two parts. The first is the canonical dimension and 
the other is the anomalous dimension. Here 

ct = + , p = 2 - ri - % (29) 

where r\ is the anomalous dimension of <p' and % is the vertex anomalous dimension. 
The former is related to the renormalization group function which is also termed the 
anomalous dimension, y(g), by 

n = Y(Sc) (30) 
where g c is the value of the coupling constant at the critical point, 

7(g) = A^ lnZ ^ (3 1 ) 

and is the wave function renormalization constant. (In (l3TT l we have temporarily 
used ji to denote the standard renormalization group scale that underlies any renor- 
malization group equation.) To determine the values of the exponents to a particular 
order in 1 /N requires solving the skeleton Schwinger-Dyson equation for the 2- 
point functions at the same order. We do not discuss that formalism here, which can 
be found in 0] [2), as our focus is rather on the evaluation of the Feynman graphs 
contributing to these equations. Though we should say that the presence of the non- 
zero anomalous dimensions in the propagators means that in 2-point functions there 
are no self energy corrections on any internal propagator as otherwise there would 
be double counting. So the number of graphs to consider is smaller than the corre- 
sponding perturbative case. 

The coupling constant at the critical point is denoted by g c and is defined as a 
nontrivial zero of the j3-function, j5(g c ) = 0. As we are working in ^-dimensions 
such a non-trivial zero exists in our theories since away from the spacetime di- 
mension where the theory is renormalizable the coupling constant becomes dimen- 
sioned. Hence the first term of the cZ-dimensional /3 -function depends on d. More- 
over, g c will depend on the parameters of the theory which in our case here is N. 
Thus g c = g c (d,N). Similarly 77 = rj(d,N) and % = %(d,N). These can all be ex- 
panded in powers of 1 /N where is large in such a way that the coefficients of 1 /N 
are of-dependent. Thence if one expresses these coefficients in powers of £ where 
d = 4 — 2e for 4 theory or d = 2 + e for the nonlinear a model, then one can 
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deduce the coefficients in the corresponding renormalization group equation to all 
orders in perturbation theory at that order in 1 /N. In this respect it is important to 
note that in the large N expansion e or e do not play the role of a regulator as they 
would do in conventional perturbation theory. 

Instead to see the origin of where a regulator is required one should consider 
the simple two loop contribution to the a self energy graph given in Figure 10 in 
coordinate space representation. 




Figure 10. Two loop self energy for <7. 

To use conformal methods one has to check the sum of the exponents at a vertex in 
coordinate space representation. From (f29} one can see that 

2a + J3 = 2 M - X • (32) 

However, from the structure of the renormalization group equation at criticality the 
anomalous dimensions 77 and % begin as (9(1 /N). More, specifically 

Thus at leading order in 1 /N the basic vertex is unique, J2). Hence at this order one 
can integrate at either of the vertices and produce the first contribution to the inte- 
gral which is V (jU — 1 , jU — 1 , 2) . The second integration is a simple chain and naively 
gives v (jJ. , jJ. ,0) . This is clearly ill-defined due to the zeroes and singularities deriv- 
ing from the F-function. However, this graph was chosen to illustrate the fact that 
the graph and indeed the theory requires a regularization in this critical point formu- 
lation. The method developed in JT] [2] was to use analytic regularization which is 
introduced by shifting the vertex anomalous dimension by an infinitesimal amount, 
A, via 

P ->• P - A . (34) 

In some respect one is in effect performing a perturbative expansion in the vertex 
anomalous dimension, 0]|2] ■ Consequently even at leading order the graph of Figure 
10 no longer has a unique vertex due to a non-zero A. Therefore, to determine the 
graph to the finite part in A requires the addition and subtraction of the graphs of 
Figure 11, Q. 
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Figure 1 1 . Subtracted graphs for computation of a two loop self energy. 

These two graphs have been chosen in such a way that their singularity structure in A 
exactly matches that of Figure 10, |2). Clearly they represent simple chain integrals 
which can be determined as 2v(a,j3 — A,2[i - a - /3 +4)v(a,ju — A, pi — a + A) 
where the singularity is clearly regularized. To complete the evaluation introduces 
another technique, which we will use later, to extract a finite term of a graph. This is 
a temporary regularization, [|2]. If one subtracts the graphs of Figure 1 1 from that of 
Figure 10, the combination is finite with respect to A which is therefore not required 
and can be set to zero. Thus one can complete the first integration at the upper vertex 
of each graph. (Without a regularization the point where one integrates in each graph 
has to be the same and thence the order of integration is important.) This produces 
v(a,a,j3) for each graph. However, each of the three subsequent chain integrals 
has a singular exponent, ji. To circumvent this the lower two propagators of all 
three graphs are temporarily regularized by a —> a — 8 where 8 is arbitrary. Thus 
the three graphs give 

[v(/i - 5,/i- 8,28) - v(a-8,/i -8,ji -a + 28) 
-v(ji-8,a-8,ii-a + 28)]v(a,a,l3) (35) 

which is clearly finite as 8 — > 0, J2). Thus to 0(A) the graph of Figure 10 evaluates 
to, @, 



2n 2 ^a 2 {a)a{^) 



j+B(li)-B(a) + 0(A) 



(36) 



r(M) 

where B(z) = y/(z) + yf(jJ. — z) for z and (jtt — z) not equal to zero or a negative 
integer and y/(z) is derivative of the logarithm of the F-function. 

A more involved example which uses many of the techniques of the previous 
section occurs in the computation of the 0(1 /N 2 ) correction to the /3 -function in 
0(N) (j> 4 theory. The relevant critical exponent is (0 which is related through the 
critical renormalization group equation to the j3 -function slope at criticality. In this 
case it has the form 

» = 2 - /* + i £ < 37 > 

and the explicit forms for co n are deduced from the part of the Schwinger-Dyson 
equations corresponding to corrections to scaling. In other words the propagators of 
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d28t are extended to 



8 ij A 



[1 



+ A\x 2 r] 



{(T(0)(T(x)) 



B 



[1 



(38) 



In principle other corrections can appear here corresponding to other exponents such 
as that for the j5 -function of the nonlinear a model but one tends to focus on one 
calculation at a time. The effect of the corrections is that to deduce co„ within the 
Schwinger-Dyson formalism all Feynman diagrams with one correction insertion on 
a propagator contribute at each particular order in 1 /N. While the 0(1 /N 2 ) expres- 
sion for co appeared in lfT31 the explicit evaluation of the contributing graphs has not 
been detailed. Thus we discuss one such diagram here as the approach can be readily 
adapted to the other graphs. It is given in Figure 12. To see that it is 0(1 /N 2 ) each 
closed loop of <p' fields contributes a factor of N and each a propagator is 0(1 /N). 
This is due to the fact that the amplitude B is 0(1 /N), IT] [2]|. As there are four of 
the former and five of the latter then this gives 0(1 /N) overall which is one factor 
of 1 /N more than the previous order graph of Figure 10. Finally, another factor of 
1 /N derives from the actual Schwinger Dyson formalism used to determine ©2- The 
double line on one <7 propagator in Figure 12 denotes the B' correction. 




Figure 12. Particular graph contributing to the (j) 4 theory /3 -function at 0(1 /N 2 ). 

The presence of such a correction means that the graph is A -finite. Moreover, since 
we only want the value as a function of d rather than d and N we can replace the 
exponents of the lines by their canonical values. If one was computing 03 then the 
anomalous dimensions of each exponent would need to be retained at 0(l/N). The 
benefit of this restriction here is that of the ten vertices eight are unique. There are 
ten integrations to do over the vertices rather than the six of the loops as we are in 
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coordinate space representation. Given this high degree of uniqueness the graph can 
be reduced rather quickly to one with fewer integrations. To do this one can use a 
variety of the rules we had earlier aside from uniqueness such as conformal trans- 
formation, unique triangle, insertion at an internal or external vertex. Ultimately one 
produces the graph of Figure 13. 



1 




Figure 13. Reduced integral of Figure 12. 

This graph cannot be reduced any further since there are no unique vertices or trian- 
gles. Though various vertices or triangles are one unit from uniqueness. Moreover, 
integration by parts cannot be used since at some point one produces an unregu- 
larized exponent, such as or jj,, or a zero in a denominator factor. In some sense 
this graph could be regarded as a master integral since it arises in several of the 
other graphs contributing to the a Schwinger-Dyson equation. Moreover, it is worth 
noting that in strictly four dimensions the propagators of the graph would all have 
unit exponents. As an aside if an interested reader has been applying the conformal 
techniques to reduce the diagram and obtains similar exponents but distributed dif- 
ferently around the diagram then it will be related to that of Figure 13 by applying 
the transformations discussed for the master 2-loop self energy. We note here that if 
a conformal transformation is applied to the graph of Figure 13 with the left internal 
point as the CL base, then that would introduce a new line from the top right internal 
vertex to the base. This illustrates comments made earlier. 

To proceed further and reduce the graph to a known function of d requires an 
integration by parts but this requires modifying the integral first. Though before this 
can be achieved safely one has to introduce a temporary regularization to handle 
hidden singularities at a later stage of the computation. This technique has been 
applied by others, J3][l6|. For our case we have choosen the regularization of Figure 
14. How one chooses the temporary regularization is not unique. However, it is 
chosen here so that after application of the integration by parts rule of Figure 7 the 
resulting four graphs have either unique vertices or triangles which are 5 -dependent 
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and which regularize any singularity after subsequent integration. For the integration 
by parts we use the top left 



[1-1-8 




li-l-S 



2/1-3 + 8 



Figure 14. Temporary regularization of previous graph to reduce it to two loop 

basic graphs. 

internal vertex of Figure 14 with the line joining the quartic vertex as the reference 
line of the rule of Figure 7. This produces the four graphs of Figures 15-18. 



[1-1-8 




H-l-S 



2ii-3 + 8 



Figure 15. First graph after integration by parts. 

All but the third have at least one unique vertex while that has a unique triangle. In 
our earlier notation the first two graphs of Figure 15 and 16 are 

v(2, 1 - 8,2/1 -3 + 8)lQi -1-8, [1-1 + 8, [1-1 + 8,^-1- 8,1) (39) 



and 



v(2,ju- 1 -8, /I- 1+5)7(1+5,2^-3,1 + 5,2^-3-5,1) 



(40) 
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As both of these are 5-finite and have no (5-singular coefficients, one can set 5 to 
zero in each. The final evaluation is by a two loop reduction formula similar to those 
of Figures 8 and 9. 



2 




Figure 16. Second graph after integration by parts. 

For the remaining two graphs of Figures 17 and 18 one has to treat them together 
due to the singular propagator exponents as will be evident. After integrating the 
respective unique triangle and vertex they combine to produce 

a 3 (l)a(n~- 5)a(2ju -3 + 5) [/(ju - l,ju- 1,1 + 5,^-1 -5,jU-l 

- I(n-l,n-l,l + 8,n- 1-25,^-1 + 5)] . (41) 



1 




Figure 17. Third graph after integration by parts. 

As the external coefficient includes a factor of 1 /5 then the quantity inside the 
square brackets needs to be evaluated to 0(8). 
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J" 




211-3 + 8 



Figure 18. Fourth graph after integration by parts. 

This is not possible exactly for both integrals. (It is for the first.) Instead since one 
only needs the 0(8) part itself one can achieve this by evaluating the integral 



/(jU- l,/x-l,l,/x-l-5,/x- 1 + 5) 



(42) 



From the two 2-loop graphs we are interested in the 0(8) term of this integral clearly 
corresponds to the piece we require. Moreover, it can be evaluated exactly using 4- 
as it then reduces to an integral to which one can apply a 2-loop recurrence relation 
similar to that of Figure 9. The final expression for the graph of Figure 12 is 



where 



a(jU - l)a 2 (2jU - 3)a(2jU -2) 
2(M-3)(M-2) 9 



h-fl 



2/i 



(M-2) 



•6/3 



(43) 



fx = yA(3-M) + V(2M-3)-y(M-l)-V(l) 
h = i/(3-M)-^(2M-3) + ^'(M-l)-'/(l) 

/ 3 = f Gu-i)-y'(i)- 



(44) 



Setting jj. = 2 reproduces the established leading order value for the wheel of 
three spokes, ifTTl . which provides a useful check. Finally, all the other contributing 
graphs are evaluated in a similar way and the full expression for ©2, after using the 
Schwinger Dyson formalism, is given in I1T31 . 



5 Future Directions 



We close the article by discussing several directions in which this approach could 
move. First, the extension of scalar field theories to non-abelian gauge theories has 
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been considered in |[T8l[T9ll20ll2Tll22l for various applications where information is 
needed on the renormalization group functions of operators in deep inelastic scatter- 
ing and the j5 -function. That approach is based on the observations of ll23ll using the 
number of quark flavours, Nf, as the expansion parameter. Rather than use the full 
QCD Lagrangian one exploits the critical point equivalence with the non-abelian 
Thirring model, 11231 . 

L = iyW + g¥fAlT?yj - i(A£) 2 (45) 

where is the covariant derivative, T" are the group generators and is the 
quark field with 1 < i < Nf. The spin-1 auxiliary field A"^ plays the role of the 
gluon in the higher dimensional theory. The triple and quartic gluon vertices of 
QCD are generated by the 3-point and 4-point functions of (l45T l with external 
legs respectively. Following the critical point analysis the propagators in a similar 
notation, but in momentum space, are 



[¥{-p)v 3 (p)) 



sab d 



p 



(46) 



where E, is the gauge parameter with the Landau gauge corresponding to E, = 1 . 
From dimensional analysis the exponents are now, |[T8l . 

a = ii + ±77 , J3 = 1 - tj - % (47) 

which means the basic vertex is one step from uniqueness. This complicates compu- 
tations in that to proceed one has to break all contributing graphs into scalar integrals 
and treat them by transformations, subtractions or use integration by parts to reduce 
them to computable cases. While it has been possible to do this in certain instances, 
OH 123122), it is not systematic. 

Since the application of the method of IT] (2) to QCD an algorithm has been de- 
veloped which allows one to exploit integration by parts. Known as the Laporta 
algorithm, QTj, it creates all integration by parts relations between integrals of a 
particular topology and then algebraically solves them in such a way that all inte- 
grals are reduced to a basis set of master integrals. Once their values are known 
by other methods then the problem is complete. In the large Nf context once one 
moves to say 0(1 /Nf) computations then graphs such as that of Figure 12 need to 
be computed in QCD. Then the solid lines would represent quarks and the springs 
would correspond to gluons. However, taking the traces over the closed loops results 
in a huge number of irreducible numerator scalar products. While the propagators 
do not have integer powers, as is the case in perturbative calculations, there appears 
to be a similarity to the problem. In other words in principle a generalization of the 
Laporta algorithm should be able to produce a reduction of the irreducible graphs to 
a set of masters. The difficulty is that the presence of non-integer propagator powers 
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means that the present Laporta algorithm would need to be modified in order to have 
a point, akin to a ground state, below which no more reductions could be possible. 
Though it is not clear under what conditions such a bottom point exists or whether 
for certain topologies or distribution of non-unit exponents it can be proved to be 
impossible. Indeed the latter point could be related to the issue of lack computability 
of a graph due to the presence of multiple-zeta values similar to (fl3T l. However, it 
seems that for the practical problem of deducing the QCD j5 -function at 0(1 /Nj) 
such an extension to the Laporta algorithm is possibly the only feasible tactic at 
present. 

Aside from this possible extension to the Laporta algorithm another interesting 
possibility is to what extent the conformal integration methods can be built into that 
algorithm to improve and speed reductions within a computer algebra programme 
for massless Feynman graphs. This may be important for higher loop topologies. 
For instance, earlier we derived recurrence relations for the two loop self energy 
topology based on the transformation deduced from the generalized uniqueness con- 
dition. While such relations are no doubt contained within integration by parts re- 
lations of the Laporta construction, that of Figure 9 is particularly useful in that 
there is no increase in the power of any propagator. Therefore, it may be possible 
to construct similar relations using conformal transformations but for higher loop 
massless topologies. Indeed such transformations are not unrelated to the symmetry 
group of the topology as has already been studied in depth for the two loop self 
energy, J6] (TJ . At the time of J6]|7| expanding a graph in terms of its group invari- 
ants was a promising approach which was complemented by later methods such as 
llllD. However, it may be worth returning to a group theory analysis for topologies 
such as that represented in Figure 13. This is because the high order expansion in 
terms of e of this and other three and four loop topologies will soon be required for 
extending QCD to^ve and possibly higher loops. In this respect another direction 
of exploration may be to study the structure of the graph polynomials of a topology. 
The transformations of |2] have been derived from a graphical approach to under- 
standing the structure of the two loop self energy graph. Understanding the effect 
such conformal transformations have on the graph polynomials of massless integrals 
may also give insight into the as yet undetermined group theory properties of higher 
order topologies. 
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